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Abstract 

The first direct measurement of electron-neutrino quasielastic and quasielastic-like scattering on 
hydrocarbon in the few-GeV region of incident neutrino energy has been carried out using the 
MINERvA detector in the NuMI beam at Fermilab. The flux-integrated differential cross sections 
in electron production angle, electron energy and Q 2 are presented. The ratio of the quasielastic, 
flux-integrated differential cross section in Q 2 for v e with that of similarly-selected ^-induced events 
from the same exposure is used to probe assumptions that underpin conventional treatments of 
charged-current v e interactions used by long-baseline neutrino oscillation experiments. The data 
are found to be consistent with lepton universality and are well-described by the predictions of the 
neutrino event generator GENIE. 
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I. INTRODUCTION 


Current and future neutrino oscillation experiments hope to measure CP violation in 
the neutrino sector by making precise measurements of z/ e (/7 e ) appearance in predominantly 
z/ ;i (z7 /t ) beams. These experiments (such as NOvAjTj, T2K[2], and DUNE (3]) consist of large 
detectors of heavy nuclei (e.g., carbon, oxygen, argon) to maximize the rate of neutrino 
interactions. They examine the energy distribution of interacting neutrinos and compare 
the observed spectrum with the predictions based on different oscillation hypotheses. Correct 
prediction of the observed energy spectrum for v e interactions requires an accurate model of 
the interaction rates, particle content, multiplicity and outgoing particle kinematics. In other 
words, there is a need for precise v e cross sections on the appropriate detector materials. 

The relatively small components of v e and V e flux in neutrino beams coupled with signif¬ 
icant backgrounds arising from the dominant interactions have led to a paucity of v e and 
v e measurements in this energy range (0.5 to a few GeV). Gargamelle[5] and T2K|[5] have 
published v e inclusive cross-section measurements at these energies, but small statistics and 
the inclusive nature of both of these measurements limit their usefulness for model compar¬ 
isons and as a basis for tuning simulations. Therefore, most simulations, such as those used 
in oscillation experiments, begin by tuning to high-precision z/ M (z/ M ) cross-section data and 
apply corrections such as those discussed in Ref. [B| to obtain a prediction for the ry(T e ) 
cross section. 

This Letter reports measurements of v e and V e charged-current quasiclastic (CCQE) in¬ 
teractions {y e n —>■ e~p and u e p —> e + n) on nucleons in a hydrocarbon target at an average 
v e energy of 3.6 GeV. Quasielastic scattering is a two-body process that is of particular 
importance in neutrino physics since it is the dominant reaction near 1 GeV, which is a 
critical energy region for accelerator-based long-baseline oscillation experiments. Though 
the incoming neutrino has an unknown energy and the final-state nucleon may not be de¬ 
tected, knowledge of the incoming neutrino direction and the outgoing lepton momentum 
vector, along with the assumption that the initial-state nucleon is at rest, are sufficient to 
constrain the kinematics. Thus the assumption that quasielastic scattering takes place on 
free, stationary nucleons is often used to extract an estimate of the neutrino energy and 
the square of the four-momentum transferred to the nucleus (E® E and Q q E , respectively). 
However, hadrons exiting the nucleus after the interaction can rcinteract and change identity 
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or eject other hadrons[7], and the complex interactions within the initial nuclear environ¬ 
ment can deform the inferred kinematics or cause multiple nucleons to be ejected by a single 
interaction [HI E]- Thus, true quasiclastic events cannot be reliably isolated experimentally. 
As an alternative, this analysis defines “CCQE-like” events to be the signal. These are 
events having a prompt electron or positron from the primary vertex plus any number of 
nucleons but devoid of any other hadrons or associated 7 -ray conversions. Both v e - and 
I7 e -induced CCQE-like events are included since the final-state e ± cannot be distinguished 
in MINERvA’s unmagnetized tracking volume. The v e have a significantly smaller flux and 
cross section relative to the zz e , though there is a small analysis selection bias favoring V e 
over v e . According to the simulation, the Xvinduced events comprise 8.9% of the selected 
sample of v e and v e interactions. In this Letter, the v e (positron) content is included when 
referring to the signal. 

The relatively high statistics in the MINERvA data set allows for flux-integrated differ¬ 
ential cross-section measurements for the v e quasiclasticlike process as well as a comparison 
of the v e and quasielastic cross sections as a function of Q q E . These measurements are 
useful for neutrino oscillation experiments seeking to quantify their understanding of the ex¬ 
pected v e energy distribution. Notably, the target medium for this analysis (hydrocarbon) is 
nearly identical to that used in NOvA and the T2K near detector, and the neutrino energy 
range of this analysis overlaps that of NOvA and DUNE. 


II. THE MINERVA EXPERIMENT 

MINERvA records interactions of neutrinos produced in the NuMI beam line |10j. In 
NuMI, a beam of 120-GeV protons strikes a graphite target and produces charged mesons 
which are focused by two magnetic horns into a 675-rn helium-filled decay pipe where most of 
the charged mesons decay producing neutrinos. For the data used in this analysis, the horns 
focused positive mesons, resulting in a beam enriched in neutrinos with a most probable 
neutrino energy of 3.1 GeV. This analysis uses data taken between March 2010 and April 
2012 with 3.49 x 10 20 protons on target (POT). 

The neutrino beam is simulated by a GEANT4-based model|ll, [12] constrained to re¬ 
produce hadron production measurements jT3b2Tj . Hadronic interactions not constrained by 
the external hadron production measurements are predicted using the Fritiof Precompound 
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(FTFP) hadron shower model [22]. The uncertainty on the prediction of the neutrino flux 
depends upon the precision in these hadron production measurements, uncertainties in the 
beam line focusing system and alignment |23j. and comparisons between different hadron pro¬ 
duction models in regions not covered by the external data. Recently, an in situ MINERvA 
measurement of purely leptonic v — e elastic scattering from atomic electrons[21] became 
available and can be used to provide a data-based constraint for the flux estimate by compar¬ 
ing the precisely predicted rate for this process with what is observed. The calculated v e + V e 
flux for the analysis in this Letter, which includes the application of the u-e constraint, is 
shown in Fig. [3] and is provided in tabular form in the ancillary material. 


xIO 9 



FIG. 1. The v e + V e flux as a function of neutrino energy from the beam simulation for the data 
used in this analysis. The u e flux is shown separately to emphasize the dominance of v e in the 
sum. 

The MINERvA detector consists of a core of scintillator strips surrounded by electro¬ 
magnetic and hadronic calorimeters on the sides and downstream end of the detector. The 
target and tracking region for this analysis is 95% CH and 5% other materials by weight. 
The triangular 3.4 x 1.7 cm 2 strips are approximately perpendicular to the beam axis and 
are arranged in hexagonal planes of three orientations, enabling stereoscopic reconstruction 
of the neutrino interaction vertex and outgoing charged tracks. The downstream electro¬ 
magnetic calorimeter (ECAL) is identical to the tracking region except for the addition of 
a 0.2-cm (0.35 radiation lengths) lead sheet in front of every two planes of scintillator. 

MINERvA is located 2 m upstream of the MINOS near detector, a magnetized iron 
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spectrometer[25]. which is used to reconstruct the momentum and charge of ^. The 
MINERvA detector’s response is simulated by a tuned GEANT4-based[IT( El] program. 
The energy scale of the detector is set by ensuring that both the photostatistics and the 
reconstructed energy deposited by momentum-analyzed beam-related muons traversing the 
detector agree in the data and simulation. The calorimetric constants used to reconstruct 
the energy of electromagnetic showers, including corrections for passive material[26] and 
algorithm-specific tuning, are determined from the simulation. Detailed descriptions of the 
MINERvA detector configuration, calibrations, and performance can be found in Refs. [56- 

M\- 

Neutrino interactions are simulated using the GENIE 2.6.2 event generator [29]. The simu¬ 
lation is used for efficiency corrections, unfolding, and background estimation. Weak interac¬ 
tion [vector minus axial-vector (V-A)] phenomenology is used for quasielastic interactions[30j 
in the simulation, with axial mass Ma— 0.99 GeV and a relativistic Fermi gas nuclear model. 
The modeled charged-current cross sections differ for v e and v jL only in the lepton mass, which 
appears in kinematic factors in the differential cross-section expressions. 


III. EVENT RECONSTRUCTION AND ANALYSIS 

Events selected for this analysis are required to originate from a 5.57-ton fiducial volume 
in the central scintillator region of MINERvA. The energy depositions in the scintillator 
strips (hits) are first grouped in time and then spatially grouped into clusters of energy 
in each scintillator plane. Clusters with energy > 1 MeV are matched among the three 
views to create tracks. The hits in each scintillator strip are recorded with 3.0-ns timing 
resolution, allowing separation of multiple interactions within a single beam spill. Candidate 
events are created from tracks whose most upstream energy deposition is in the fiducial 
volume and which do not exit the back of the detector, as such highly penetrating tracks are 
overwhelmingly muons. All tracks passing the criteria above are tested as e ± candidates. 
Hits are considered if they fall within a region that consists of the union of two volumes: 
a cylinder of radius 50 mm extending from the event vertex along the track direction and 
a 7.5° cone with an apex at the event vertex (origin of track) and a symmetry axis along 
the track direction. Hits are associated with the cone as it extends through the scintillator 
tracker and ECAL; the collection of hits ceases when a gap of three radiation lengths is 
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encountered that is devoid of hits. The hits in this cone “object” are examined using a 
multivariate particle identification (PID) algorithm. This technique combines details of the 
energy deposition pattern both longitudinally (mean dE/dx and the fraction of energy at 
the downstream end of cone) and transverse to the axis of the cone (mean shower width) 
using a fc-nearest-neighbors algorithm|31j. For those candidate events deemed consistent 
with an electromagnetic cascade, electrons and positrons are separated from photons by 
demanding the energy deposition near the upstream end of the cone be consistent with a 
single track rather than the two particles expected from photon conversion to e + e _ . The 
discriminant used for this separation is the minimum energy in a sliding 100-nmi window 
along the axis of the cone, in 20-mm steps, from the event origin up to 500 mm (about 
1.2 radiation lengths). This technique reduces the possibility of bias introduced by nuclear 
activity near the interaction point |32j. Cone objects surviving to this point are considered 
to be electron (or positron) candidates. 

The next stage of the analysis requires the topology of the event to be consistent with 
v e CCQE-like. Events containing tracks consistent with charged pions or muons or events 
with electromagnetic activity outside of the electron candidate cone object (such as might 
be expected in the presence of a 7r° decay) are removed by a cut on the “extra energy ratio” 
variable T. This quantity represents the relative amount of energy outside the electron can¬ 
didate cone to that inside the electron candidate cone. Hits within a sphere of 30-crn radius 
about the interaction vertex are ignored when calculating T to reduce the contribution from 
low-energy nucleons which are potentially not well simulated [32]. Events at large T are 
removed from the sample. The cut in T is a function of the total visible energy of the event 
and was tuned using simulated events. In addition to the T cut, Michel electron candidates 
from the tt —> /i —> e decay chain are rejected via timing and their spatial proximity to track 
ends. 

Finally, events are retained in the sample only if they have a reconstructed electron en¬ 
ergy E e greater than 0.5 GeV and a reconstructed neutrino energy E® e less than 10 GeV. 
The lower bound excludes a region where the expected flux of v e and z/ e is small and the 
backgrounds are high. The upper bound eliminates events in the region of large flux uncer¬ 
tainty. 

The reconstructed electron energy distribution of the 2105 selected v e CCQE-like candi¬ 
dates is shown in Fig. [2] for both the data and the simulated event samples. The simulated 



FIG. 2. The reconstructed electron energy distribution after all selection cuts and after constraining 
the backgrounds using sidebands in the data. The errors shown on the data are statistical only. 

sample is broken down by process according to the GENIE event generator and is 52% 
pure signal events. The primary source of background in the selected sample arises from 
i/ e -induced non-CCQE-like events. The second largest background comes from incoherent 
neutral current (NC) 7T° production. Coherent 7T° production and neutrino-electron elastic 
scattering also contribute to the final sample. 

The sizes of the backgrounds in Fig. [2] are constrained by two sideband samples. The Erst 
sideband consists of events at larger T, which is enriched in inelastic backgrounds from v e 
interactions and incoherent events containing 7r°. The other sideband, dominated by v e CC 
inelastic events, consists of events with Michel electron candidates (where the Michel electron 
was typically produced via the decay chain of a charged pion). The normalizations of the 
u e inelastic and incoherent 7r° backgrounds are varied in order to find the best overall fit of 
simulation to the data in the reconstructed electron angle and reconstructed electron energy 
distributions in each sideband sample. Since, according to the simulation, the sideband in T 
contains some signal events, the procedure is iterative. The background scale fit is done and 
the signal is extracted and used as a constraint for a new background scale fit. This is done 
until the background scale factors stabilize (two iterations). After this procedure, the fitted 
scale factor for the normalization for the u e inelastic category is found to be 0.89 ± 0.08, 
while that for the incoherent 7T° processes is 1.06 ± 0.12. The neutral-current coherent 
pion production is scaled down by a factor of 2 for pions with energies below 450 MeV in 
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the simulation to bring the GENIE charged-current coherent charged pion production into 
agreement with a recent MINERvA measurement [33]. Subsequent to these constraints, the 
scaled backgrounds in the signal region are subtracted from the data. 

An excess of photon or 7r°-like events in the data relative to the simulation was observed 
in the distribution of energy deposited in the upstream part of the electron candidate cone, 
as characterized and described in detail in another paper{31]. Models of single photon or 7r° 
production consistent with the observed excess were evaluated and found to have little effect 
on the background in the signal region of this analysis. Nevertheless, a 7r° background fitted 
to the excess is added into the simulation and contributes (negligibly) to the background 
subtraction. 

The flux-integrated differential cross sections in electron energy E e , angle 9 e , and four- 
momentum transfer Q q E are calculated in bins i as a function of sample variable £, with e 
representing signal acceptance, <f> the flux integrated over the energy range of the measure¬ 
ment (or over the bin i, in the case of the total cross section), T n the number of targets 
(nucleons) in the fiducial region, Aj the width of bin i, and U t j a matrix, derived from the 
simulation, correcting for detector smearing between bins i and j in the variable of interest: 


/ dcr\ 1 

V di Jei$T n Ai 


- x Uij (V data - A bknd pred ). 


( 1 ) 


E^ E and Q q E are calculated from the lepton kinematics alone using the approximation of 
a stationary target nucleon. Unfolding to correct for detector effects in the four variables is 
done using a Bayesian technique [35] with a single iteration. 

The systematic errors considered arise from the primary neutrino interaction model, the 
flux model, and the detector response to particle activity. The errors on the flux are deter¬ 
mined as discussed earlier. At the focusing peak, i.e., those neutrinos most relevant for this 
analysis, the v e flux arises from muons from pion decays. The errors in the primary neu¬ 
trino interaction model are evaluated via the reweighting of events by varying the underlying 
model tuning parameters according to their uncertainties. The parameters varied in this way 
include the shape and normalization for elastic and resonance productions, nuclear model 
parameters principally affecting the deep inelastic scattering, and parameters which con¬ 
trol the strength and behavior of the final-state interactions. Contributions to the detector 
response systematic error were determined by varying the energy scale for electromagnetic 


10 





FIG. 3. Flux-integrated differential v e CCQE-like cross section versus electron energy (left) and 
electron angle (right). Inner errors are statistical; outer are statistical added in quadrature with 
systematic. The band represents the statistical error for the Monte Carlo curve. 

interactions, the parameter used in Birks’ law, the photomultiplier tube cross-talk fraction, 
the Michel electron reconstruction energy scale, and the detector mass. The largest system¬ 
atic errors contributing to the cross-section results presented here are due to the detector 
response, the interaction model, and the flux model, with each contributing a fractional 
uncertainty of less than 10%. The overall systematic errors are typically in the 10%-15% 
range, which is sufficiently small for the results presented here to be statistically limited. 

The flux-integrated differential v e CCQE-like cross sections versus electron energy and 
angle are given in Fig. [3j for both the data and the POT-normalized Monte Carlo samples. 
The analogous distribution in Qq E is given on the left side of Fig. |4j The measured cross 
sections and covariances are provided in tabular form in the ancillary material. The sim¬ 
ulation appears to underestimate the width of the electron production angle and exhibit 
a harder spectrum in Qq E . However, these differences are not significant when correlated 
errors, such as the electromagnetic energy scale, are taken into account. 

In order to compare directly the measured differential cross section for v e and zy, inter¬ 
actions on carbon as a function of Q q E , an analysis similar to that described in this Letter 
was performed in terms of a CCQE signal (rather than CCQE-like), as specified by the 
GENIE event generator, which can be compared directly to previously published MINERvA 
results[32], The selection cuts for the u e events were adjusted slightly to ensure the energy 
range of included events agreed with that of the analysis. The ratio of these two results 
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FIG. 4. The flux-integrated differential u e CCQE-like cross section versus Q q E (left). Inner errors 
are statistical; outer are statistical added in quadrature with systematic. On the right is shown the 
ratio of the MINERvA v e CCQE differential cross section as a function of Q q E to the analogous 
result from MINERvA for v,, [32|. In both figures, the band represents the statistical error for the 
Monte Carlo curve. 


and the corresponding ratio of the Monte Carlo predictions are given on the right in Fig. |4} 
The data for the differential cross section for v e CCQE interactions agree within errors with 
that for v jt CCQE interactions. (Some of the uncertainties evaluated in this analysis, such 
as the electromagnetic energy scale, result in Q 2 -dependent changes to the data distribution 
shape. These can cause trends similar to the upward slope in Fig. [4} When accounting for 
these correlations, the shape of the data curve is consistent with the shape of the GENIE 
prediction within la.) 


IV. CONCLUSIONS 

This Letter presents the first exclusive measurement of the flux-integrated differential 
cross section for v e CCQE-like interactions and thus provides the first data for directly test¬ 
ing and tuning models of a critical channel for accelerator-based oscillation experiments. 
The flux-integrated differential distributions of the cross section in E e , 9 e and Q q E agree 
with the expectation from lepton universality. A direct comparison, in the same detector, of 
the differential flux-integrated cross section of v e CCQE interactions to that for CCQE 
interactions as a function of Q q E also shows good agreement. Collectively, these measure- 
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ments constitute an important first test of the common assumption made by oscillation 
experiments that z/ M cross-section data can be applied to models of v e CCQE interactions. 
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